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A perturbative expansion of Loop Quantum Cosmological transitions amplitudes 
of Bianchi I models is performed. Folio-wing the procedure outlined in [1, 2] for 
isotropic models, it is sho-wn that the resulting expansion can be -written in the form 
of a series of amplitudes each -with a fixed number of transitions mimicking a spin 
foam expansion. This analogy is more complete than in the isotropic case, since 
there are no-w the additional anisotropic degrees of freedom -which play the role of 
'colouring' of the spin foams. Furthermore, the isotropic expansion is recovered by 
integrating out the anisotropies. 
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^; I. INTRODUCTION 

i^ . An important open problem in Loop Quantum Gravity (LQG) [3-6] is to obtain a -well 

defined method of perturbatively computing its dynamics. The covariant approach given by 
_ Spin Foam Models (SFM) [4, 7, 8] provides an avenue to obtain such a method, ho-wever 

^^ • there are still several open issues regarding its precise relation to the Hamiltonian theory. 

CN . If SFM and canonical LQG are to be the covariant and canonical descriptions of a single 

quantum theory of gravity, one should be able to derive one from the other. We are still 
_ far from rigorously establishing such connection, ho-wever important progress has arisen in 

Q ' recent years. These include, in the covariant — ^ canonical direction, the derivation of the LQG 

O . Hilbert space as -well as the spectra of geometrical operators [9-11], and, in the canonical 

— )■ covariant direction, the extension of the EPRL amplitude to arbitrary valent spin foams 
vertices thereby allowing general histories of graphs [12]. The latter direction leads also 

^ ■ to the picture of regarding spin foams as spin net-works histories [13]. This interpretation 

however, has not gone beyond the heuristic level (as far the full four dimensional theory is 
concerned; in the three dimensional case the connection between the canonical and covariant 
descriptions is well established [14]). 

Recently, the canonical — )• covariant direction was analyzed at the symmetry reduced level 
of homogeneous and isotropic cosmology [1 , 2] . Here we extend that analysis and consider 
a non-isotropic cosmological model. The additional degrees of freedom allow for a richer 
discussion than in the isotropic case. In particular, besides the 'vertex expansion' present 
in the Friedman-Robertson- Walker (FRW) case, there are additional sums over the extra 
parameters, which are interpreted as giving a 'colouring' of the graph, thus strengthening 
the analogy with spin foams. 
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The aim here is to obtain a 'sum over histories' description within Loop Quantum Cosmol- 
ogy (LQC) [15]. Our construction is then different from that of 'spinfoam cosmology' [16, 17], 
where cosmological amplitudes are obtained using SFM as the starting point. An interesting 
question which we do not address here is whether there is any precise relation among the 
two constructions. Let us also mention that the idea of looking for a canonical-covariant 
connection at the homogeneous level has appeared before in the context of Plebanksi the- 
ory [18]; such approach could shed light into the previous question of relating 'cosmological 
spinfoams' with 'spinfoam cosmologies'. 

The paper is organized as follows. In Section II we introduce the model we will work 
with, namely the quantum Bianchi I cosmology with a massless scalar field, as obtained 
by Ashtekar and Wilson- Ewing in [19]. In Section III we construct the sum over histories 
description of the model. We outline how individual amplitudes are to be calculated and 
illustrate the procedure for a simple history. In Section IV the vertex expansion of the 
FRW [1] model is recovered by integrating out the anisotropics of our vertex expansion. We 
finish the paper with a discussion in Section V. 



II. LOOP QUANTUM COSMOLOGY OF THE BIANCHI I MODEL 

We are interested in the Bianchi I cosmological model, which is the simplest non-isotropic 
homogeneous cosmology, coupled with a massless scalar field 0. As in the isotropic case, 
one can fix a fiducial 3-metric ds^ and choose Cartesian coordinates on the spacial slice such 
that dsp = (^z1 + d2;| + dz^. The physical 3-metric is then determined by three directional 



scales factors, ai,a2 and as. 



ds^ = a\ Azl + al dzj + al dz^ . (2.1) 



The Hamiltonian analysis requires one to choose a fiducial cell V, for which we take the 
rectangular prism < Zi < Li, for each direction i = 1,2,3. The physical volume of the 
cell is then given hj V = \a1a2as\L1L2L3. Note that the choice of the fiducial cell V (i.e. 
the choice of Li, L2 and L3) is arbitrary, and one has to ensure that physical results are 
insensitive to that choice (see [19] for further discussion). 

When one goes to the quantum theory [19], it is convenient to work with a new set of 
variables, (Ai,A2,f), defined by 



X _ sgn{ai)^\a2a3\L2Ls 



_ sgn{a2)^/\a3al\L3Ll 
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2A1A2A3. (2.4) 



Here io is the square root of the 'area gap' A = 4\/37r7 £pj, 7 is the Barbero-Immirizi param- 
eter, and A3 is defined in a similar way as Ai and A2. In this representation, the gravitational 
Hilbert space "Hf^^^ consists of functions ^(Ai, A2,f) with support on a countable number 
of points and with finite norm ||\l'|p := XIa a- v \'^{-^1j-^2,v)\'^ < 00. The matter Hilbert 
space is the standard one: "HJ^^jJ** = L^(M, d0). The total kinematical Hilbert space is thus a 



tensor product Tikm = 'HkiT ® ^Sn" ^^^^ ^^ usual in LQC, the dynamics of the system are 
encoded in the constraint equation 

-C^ = dl^ + e^ = , (2.5) 

where O is a symmetric operator that acts on the gravitational part. As in [19], we restrict 
attention to the 'positive octant' {v, Ai, A2 > 0). The action of takes the form, 

(6^) (Ai, A2, v) ~V^[{v + 2)v/^7T4 ^+(Ai, A2, v)-{v + 2)^^+{X,, A2, v) 

-{v- 2)y^^o-(Ai, A2, v) + {v- 2)^/\v - 4| ^4 (Ai, A2, v)'\ , (2.6) 

where \E'^4 are defined as: 

/f ±4 f ±2 \ /f±4 

^^(Ai,A2,t;)=^ ^— -Ai, A2,t;±4 +^ ^— ■ Ai, A2,i; ±4 

\t> ± 2 V J \v ± 1 

^ fv±2 ^ v±A ^ , A T /^^±2 ^ ^ 

+ ^ Ai,^— ■A2,t;±4 +^ Ai,A2,t;±4 

\ V V ±2 J \ V 

/ V ±2 \ / V ±A \ 

+ ^ ( Ai, -^^ ■ A2, t; ± 4 j + ^ ( Ai, ^^^ ■ A2, t; ± 4 j , (2.7) 

and 

^^ Ai, A2, v) = -^[ Ai, —— ■ A2, t; + ^ Ai, A2, V 

\ V V ±2 J \ V 

+ '^(^--Xi,- A2,t;) +^(^--Ai,A2,t; 

\t> ± 2 V J \v ± 2 

+ ^ (ai, ^^ ■ A2, t; j + ^ f Ai, ^^ ■ A2, v\ . (2.8) 

As noted in [19], the operator preserves the subspaces "H^ C "H^^^^ of states whose support 
lies on the lattice 

v = e + AZ. (2.9) 

These superselection sectors have the same form as in the isotropic case [20], and, as done 
there, we will restrict to the physically interesting e = sector (the one that contains the 
classically singular value v = 0). The space we finally work with is the space of vectors with 
support on the positive octant and the e = lattice, which we denote by Hq^^ . 

We now introduce a different representation of the space 'Hq'^'^ C "Hf^^^, by changing 
coordinates (Ai,A2,w) — )■ (n, x) where 

n ■= -V eN , (2.10) 

x=(xi,X2) := (logAi,logA2)GM'. (2.11) 

In this representation, states are described by functions \E'(n, x), which again have support 
on a countable number of points and have a finite norm ||\E'|p = Yly:v l^l^?"^)!^ < '^^- 



They represent the components of the state in a {|n, x)} basis, which is characterized by 
eigenvalues of Ai, A2 and V = 4:n'yip/o^i^2^3 as follows 

V\n,^) = 87r74/on|n, x) , (2.12) 

Aj|n,x) = logXj|n,x), i = l,2, (2-13) 

with their normalization given by the product of Kronecker deltas: 

{n, X |n, x) = (5„,„/(5x,x'- (2.14) 

In the subsequent sections, we will regard "Hj]"'"'' as a tensor product of the 'volume' factor 
and the 'anisotropy' factor: 

n+++ = -Hv^nx, (2.15) 

|n,x) = \n) (g) |x) . (2.16) 

Within this splitting, O is expressed as a sum of the tensor product of operators acting on 
T-Lv and "Ha, 

6 = ^|n + l)(n| (8 ©(„+!)„ + \n){n\ (g) ©„„ + \n - l){n\ (g) 0(„_i)„ . (2.17) 

n 

The form of the operators acting on the anisotropy factor is quite simple: they are composed 
of translations in the x plane of lengths 

, , 2n ± 1 , , 

a^=log^— , (2.18) 

which depend on the volume n. If we write the operator generating the translation by a„ in 
the Xi direction as 

(e-nPi^) (n, xi, X2) := ^(n, x^ + a„, x^) , (2.19) 

and similarly for translations in the X2 direction, the operators acting on Tix take the form 

Qnn = 27rG' [n{2n + 1) [cosa+pi + cosa+p2 + cos(a+j92 - a'^pi)] + n{2n - l)[a'^ -)■ a^^]) , 

(2.20) 
e„±i „ = -7TG^n{n±l){2n ± 1) fe~'''"P' + e'""±V2Pi ^ ^-iaipi^i<±r/2P^ + Pi ^ P2) ■ 

(2.21) 
One can easily verify that the action of Eq. (2.17), when written in terms of the original 
representation, reproduces Eq. (2.6). 

Let us conclude this section by mentioning a remarkable property of the operator. As 
found in [19], one can recover the FRW cosmology by 'integrating out' the anisotropics of 
the Bianchi I model. Specifically, it was shown that there is a projection map from the 
Bianchi I space to the FRW space defined by 

J]^(n,x;0) = M/^^W(n;0), (2.22) 



in which the O operator is mapped to the O^^^ operator of the FRW model^ namely, 

J2 6>*(«, x; 0) = e''^w^^^W(n; 0). (2.23) 

X 

We will later see how this projections holds order by order in the vertex expansion. 

III. SUM OVER HISTORIES 

As in [1, 2], the natural object on which to construct a sum over histories expansion is the 
physical inner product between 'initial' |[nj,Xj,0j]) and 'final' |[n/,x/,0/]) physical states. 
This inner product is constructed from a group averaging formula involving the kinematical 

states |nj,Xj,0j) and |n/,xj,0j), 

([nj,xj,0j],[ni,Xi,0j) = 2(nj,xj,</)/|/^^dae^"^ Ip^l |ni,Xi, 0^), (3.1) 

(the \p^\ term is there so that the normalization agrees with the one used in [2]). As in the 
FRW case [1, 2], a key simplification comes from the fact that the constraint C is a sum of 
two commuting pieces that act separately on "HJJJn" and 'H^^^ ■ Consequently, the integrand 
of Eq. (3.1) splits into a matter and gravitational factors: 

2(n/,x^,0^|e^"^ |p<^||n„x„0,) =2(0^|e*"^^b<^||0,)(n^,x^|e-^"®|n„x,). (3.2) 

The matter part can be easily evaluated as, 

2(0/|e^""' b0||0.) = 2/|fe^"P'e^P^(^/-<^>)|p^| . (3.3) 

The non-triviality of Eq. (3.1) lies in the gravitational part, (n/,x/|e^*"®|nj,Xj). Following 
the strategy depicted in [1], we will express such term as a sum over histories. This can be 
achieved by observing that the term has the form of a matrix element of a fictitious evolution 
operator e~*"®, with playing the role of Hamiltonian and a that of time. 

Once the gravitational factor is written as a sum over histories, the idea is to perform 
the integral over a for each history separately, obtaining at the end a sum over histories 
expansion of the physical inner product. These steps will be discusses in the following 
subsections. 

A. Sum over histories for the gravitational amplitude 

To construct a 'sum over histories' expansion of the gravitational amplitude 
(uf, x/|e^*"® |nj, Xj), one would proceed with a Feynman-like procedure of dividing the 'time' 
a into A^ steps of length e = a/N, inserting a complete basis in between each factor, and 
finally taking the A^ — ?■ oo limit. In [2] it was shown (in the FRW context, but the result 
is generic for any discrete labeled basis) that the resulting limit is equivalent to a specific 



^ See [21] for an alternative projection that produces isotropic states, but not the 0^^*-^ associated with 
the j^-quantization procedure. 



perturbative expansion of the 'evolution' operator under study. We will use this result here 
to construct the sum over histories directly from the perturbation series. 

The starting point in such a derivation is to write the fictitious Hamiltonian as an 
'unperturbed part' Bq plus a 'perturbation' 61, 



9 = 00 + 01 . 



(3.4) 



In a spin network/spin foam picture, the above splitting would correspond to a graph 
preserving piece, 0o, plus the remaining graph changing part, 0i. In our case, we choose to 
interpret the label n as containing the information of the 'graph' , and the remaining x label 
as the colouring of the graph. Thus, 0o and 0i are respectively diagonal and off-diagonal 
in n. In the tensorial notation used in Eq. (2.17), these operators are given by. 



©0 = ^\'n){n\ ®Qnn , 
n 

©1 = ^ |n + 1) (n| (g) 0(„+i)„ + \n-l) (n| ® 0(„_i)„ . 



(3.5) 
(3.6) 



The construction now follows as in the FRW case [2], where the same label was used to 
trigger the transitions. Using standard perturbation theory in the interaction picture, the 
transition amplitude is be written as 



(n/,x/| e 



'ictOl 



|ni,Xi) = (n/,x/| 



M=0 



\M 



drju . . . / dri 
^0 



-i(a-rM)So0 g-«(rM-T"M-i)0()0 p~i{T2-ri)@a r\ -iriQo 



(3.7) 



The M-th term of the sum generates all histories with M transitions. These histories are 
obtained by inserting M— 1 identities in the form 1 = ^^ (|nm)(^m|)'8lA, m = 1, . . . , M— 1 
next to each 0i factor. This results in a sum over a sequence of volumes (no,ni, . . . ^um) 
(with tiq = rii and Hm = Uf held fixed), given by 



where 



A{nM, 



{nf,yif\ e *"® \n. 



;,Xi) = ^^ ^ (xj|A(nM,...,no;a)|xi) 

M=0 riM-l, •••)"! 



(3.^ 



(3.9) 



A(nM,---,no;a) := (-i) / ^tm ■ ■ ■ I dri A(nM, • • • ,no; tm, • • • , n; a) , 

Jo Jo 

,no; tm, . . . ,ri; a) := Anf^{a - rM)KMnM-i^nM-i(^M - tm-i) ■ ■ ■ Vnino-4„o(ri) , 

(3.10) 



with An (t) and Vn'n defined as, 



A„(r) :=e-*"«"" (3.11) 

K..„ :={®""' '^tl (3-12) 

Note that all the factors in Eqs. (3.9 - 3.12) are operators on Tix , whilst the actual 
amplitude, Eq. (3.8), involves matrix elements of the operator defined by Eq. (3.10). Note 
also that the only sequences entering in the sum are such that rim = "^m-i ± 1, -m = 1, . . . , M. 
In particular, for M fixed, there is a finite number of terms. 

The construction above has the same form as in the FRW case. The distinction however 
lies on the fact that there are additional degrees of freedom, given by the anisotropies x. 
However, in the description given so far, intermediate anisotropies do not appear since they 
are implicitly 'summed over'. To make these additional sums explicit, we insert identities 
in the form 1a = J^x^n l^m)(xm| and 1^ = J2y^ |ym)(ym| to the right and left of the 
Arimijrn+i — Tm) Operators in Eq. (3.10). The gravitational amplitude (3.7) then takes the 
form, 



(n/,x/| e *"® |ni,Xj) = 

oo 

^5]] 5Z ^(yM,nM,XM,yM-i,nM-i,---,yo,no,xo;a) (3.13) 

yo,---iyM-i 

where now we have, on top of the 'graph history' sum (given by the volume sequence), a 
sum over all possible 'colourings' for each 'graph history'. The amplitude for such history is 
given by 

A{yM, um, xm, Ym^i^um-i, • • • , yo, ^o, xq; a) := (3.14) 

pa i>T2 

{-i)^ / dTM . . . / dri A(yM, Um, xm, yM-i, ^m-i, • • • , yo, ^o, xq; Tm,---, n; a) , 



where 



Mym, nM, Xm, yM-i, nM~i, • • • , yo, ^o, Xq; Tm, ■ ■ ■ ,Ti;a) := (3.15) 



and 



A„yx(r) := (y|A„(r)|x) (3.16) 

K'xny := (x|K'n|y) • (3.17) 

Note that, in spite of their appearance, the sums over intermediate anisotropies in (3.13) 
are well defined. Although in principle the anisotropies can take any real value, in practice 
only a countable subset of the real numbers is involved in the sum (the amplitude vanishes 
elsewhere). More details on this point are given in Section IIIC 



Xi X2 



FIG. 1: The Bianchi I model can be written in terms of three parameters (n, xi, X2)- n dictates the 
volume, whilst xi and X2 are analogous to the spin labels of LQG and are represented here as two 
edges. Note that this is not intended to be a true LQG graph, rather it is a pictorial description 
of the degrees of freedom that our states have. 

B. Vertex expansion of the physical inner product 

We now use the above construction to obtain an expansion for the physical inner product, 
Eq. (3.1). First, one rewrites the integrand of Eq. (3.1) as in Eq. (3.2) and then the gravita- 
tional factor is written using the expansion given in Eq. (3.13). One then interchanges the 
integral over a with the sums over M and the intermediate labels, to arrive at a 'sum over 
histories' expansion of the physical inner product, 

oo 

([n/,x/,0/],[ni,Xj,0i]) = J^ J^ J^ A(yM,nM,XM,yM-i,nM-i,---,yo,no,xo) 

M=OnM_i,...,rai xi,...,xm 
yOv:yM-l 

oo 

= 5Z 5Z ^(«M,---,^o;Xi,x/;0i,0/) (3.18) 

where 



AiYM, Um, Xm, Ym-i, riM-i, • • • , Yo, no, Xq) := 2 J^ e^P0(^/-^.) \p^\J^^da (3.19) 

2 

e*"P.*yl(yM, um, xm, Ym-i, ^m-i, • • • , Yo, no, xq; a) . 

Pictorially, we can represent the expansion as follows. First, we represent the gravitational 
ket |n, x) as depicted in Fig. 1. A 'history' with one transition in n — )■ n is then represented 
in Fig. 2. 

Note that here the analogue to spin foams is not exact, since in a spin foam, the spin 
labels on each face of the triangulation are constant. In the expansion derived here there is 
non-trivial dynamics for the 'spins', Xj, even with a fixed 'graph' n. In general there are two 
distinct labels for a face, those at the beginning, x, and those at the end, yi, and Yi 7^ x,. 
Also, in a spin foam there is no restriction on the spins across a vertex, whereas in our case 
Yi and Xj_|_i (the 'spin' labels on either side of the {i + l)*'^ vertex ) are very closely related, 
by the form of Eq. (2.21). Thus, although the analogue is not complete, the form of the 
expansions are qualitatively the same. 




FIG. 2: The history of the spin network state i.e. the spin foam analogue of our vertex expansion, 
with M = 1 is represented pictorially. The initial state has volume n and anisotropies Xj. The 
state then 'evolves', keeping n constant but allowing the anisotropies to vary. Eventually there 
is a transition to a new volume, n. The anisotropy before and after the transition are yi and Xj 
respectively. The final state has labels n and y^. The dashed lines indicate that the 'spin' labels (xj) 
evolve along each of the constant volume pieces. The amplitude for such process, in the notation 
of Eq. (3.19), is given by yl(y,n,x,y, n,x). 

C. Histories amplitudes 



We now discuss how path amphtudes appearing in Eq. (3.18) can be calculated and 
illustrate the procedure in the simplest case. The first step is to evaluate the matrix element 
{Kf\A{nM, ■ ■ ■ ,no;TM, ■ ■ ■ ,Ti;a)\:s.i) of the operator defined in Eq. (3.10). This operator 
consist of compositions of operators A„(r) and Vn'n, given respectively in Eq. (3.11) and 
Eq. (3.12), which themselves are constructed from translations in the x plane. Because of 
the translation invariance, we will consider the matrix elements between (x| and |0); the 
original matrix element is then recovered by the substitution x — )■ xj — Xj. 

Let us discuss the structure of the operators in more detail. As already noted, Vn'n 
vanishes unless n' = n ± 1, in which case it is given by Eq. (2.21). It consist of an overall 
factor times six simple shifts involving lengths of a^ and 0'n±i/2- "^^^ ^ni^) operator is the 
exponentiation of (— ir times) the operator given in Eq. (2.20). It can be factored into two 
terms, 

(3.20) 



where 



A 



r 



A„(r) = A+(r)A„(r) , 

- jr27rG n(2?iit 1) cos a„ pi +cos a„ p2+cos{a„ P2^in Pi ) 



(3.21) 



which only involves shifts with step-size a^. The total operator A(nM 



is then a product of operators involving shifts of lengths of with j 



■■■,no;TM,---,Ti;a) 
-J ...... ^ n^ or j = n„± 1/2 

{m = 0, . . . , M). As a result the matrix element between (x| and |0) vanishes unless x lies 
in the 'lattice' generated by the of steps. The computation simplifies by selecting among 



the flj s, a set of independent (i.e. incommensurate) lengths that generate the 'lattice'. 

Let us illustrate the situation by considering the simplest path, namely the M = case. 
In this case the amplitude is given by the matrix element (x|A„(r)|0). As already noticed, 



10 

this operator contains two step-sizes, a^ and a~, and so the nontrivial matrix elements occur 
whenever x hes in the lattice ^, 

X = k+a+ + ICa- ; k+ = (A;+, A;+) G Z^, k" = {k^^k:^) E 1? . (3.22) 

From the definition of a^ and a", Eq. (2.18), one can check that these two numbers are 
always incommensurate, and so they form an independent set of generators of the lattice. 
That is to say, a point in the lattice defined by Eq. (3.22) is uniquely decomposed into its 
a^ and a~ components, i.e., if k+a^ + k~a„; = l^a^ + l^a^ then k^*" = l^*". 

Thus, the kets |x = k"'"a^ + k~a^), form a basis of the subspace of vectors which give 
a non- vanishing matrix element. This space has the structure of a tensor product of two 
copies of L^(Z^), 

|x = k+a+ + k-a;) = |k+)® |k-) (3.23) 

where |k+) and |k^) are viewed as basis of two abstract L^(Z^) spaces. 

Viewed in this way, the A^(t) operators act separately on each L^(Z^) factor: 

(x = kX + k-«;l^n(r)|0) = (k+|A+(r)|0)(k-|A;(r)|0) . (3.24) 

Each factor is the matrix element of a (translation invariant) operator in a single L^(Z^) 
space, and so can be evaluated by taking the Fourier transform. Using Eq. (2.20) one finds 



j^iMi/^XlQS / ^ / ^gife±ei+ife±92g-i27rGn(2n±l)r(cosei+cos6»2+cos(6»2-ei)) (3 25) 



27r io 27r 

In the general case of a path with M transitions and g independent generators, the space 
of vectors giving a non-vanishing matrix elements will have now the structure of a tensor 
product of g copies of L^(Z^). Because the vertex Vn'ri is the sum of six terms, one will 
generically have a total of ~ 6^ terms, each of them involving g Fourier integrals of the 
type described above. The expressions for these integrals can be directly read off from 
Eqs. (2.20) and (2.21). 

After the evaluation of the matrix element (xj|A(nM, • • • , "^o! ^m, • • • , Ti; a)|xj), one has to 
perform the r integrals in Eq. (3.10), and the p^ and a integrals in Eq. (3.19). In our M = 
case, there are no r integrals to perform, and the a integral can be done if interchanged 
with the Fourier integrals of Eq. (3.25). This gives a Dirac delta, which in turn allows one 
to evaluate the integral over p^. The result is, 

54/ 



A(n;xi,x/-;0j,0f) = / gieifci++e2fcj-+e3fcr+e4fc-7g»(</'f-<^»)V2^i^^^/(2«+i)fc(ei.e2)+(2«-i)fe(e3.e4) ^ 

Jd {^t^Y 

(3.26) 

where xj — Xj = k+a^ + k^a^, h(Q\^ 62) = cos^i + cos ^2 + cos(^2 — ^1) and the domain of 



^ Notice that the lattices we will referring to are not the usual ones (as for instance a square lattice a7? C K^) 
where points form a grid. Rather they fill in the entire plane. For instance one can show that the points 
of the lattice defined by Eq. (3.22) form a dense subset of M^ (because a^ and a^ are incommensurate 
numbers). 
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integration is 

^ = {(^1, ^2, ^3, 0^) e [0, 27r)4 / {2n + l)/i(^i, ^2) + (2n - l)/i(^3, ^4) > 0}. (3.27) 

So far we have discussed histories amphtudes of the form A(njv/, • • • ,'^o;Xj,xj; 0j, 0j), 
where intermediate anisotropies are aheady summed over. Let us now discuss briefly 
how one would compute the amplitude for an individual 'colouring' of such a history, 
^(Ym, ^m^ xju, yM-ii^M-i^ ■ ■ ■ , Jo, "^0) xo). The building blocks in this case are the ampli- 
tudes A„yx(r) and \4'xny, Eq. (3.16) and Eq. (3.17). The first amplitude coincides with the 
M = case discussed above. We thus have that A„yx(r) vanishes unless y— x = k"'"a^+k~a~ 
in which case it is given by Eqs. (3.24) and (3.25). On the other hand, the value for l^'xny 
can be easily read off from its definition: it vanishes unless n' = n ± 1 and y — x lies in one 
of the following six points, {(a^, 0), (a ^±i/2, 0), ( «„. ^^±1/2)' (0' ^n)' (0' ^^±1/2)' («n±i/2' «« )}, 
in which case it takes the value —iiG^n{n ± l)(2n± 1). After multiplying by the remaining 
factors in Eq. (3.15), one has to perform the same integrals as previously, in this case given 
in Eq. (3.14) and Eq. (3.19). 

D. Vacuum case 

The presence of matter entered only at the very end of the construction. If we did 
not have matter at all, we could still follow the same procedure and arrive at the vacuum 
equivalent of Eq. (3.18), 

00 
([n;,x;],[n„x,]) = Y.[ E ^™'""'"(nM,...,no;x„x;)] . (3.28) 

M=0 nM-l,---i"l 

The difference between the matter and vacuum cases lies only in the form of the amplitudes, 
which in the vacuum case are formally given by, 

A^--"-(no, . . . ,nM;x„x;) = /^^d« (x^|A(no, . . . ,nM\ «)|x,) . (3.29) 

These amplitudes can then be evaluated following the strategy given in the previous section, 
the only difference being the absence of the final integral over p^. It is not obvious whether 
the integral in Eq. (3.29) converges for all paths thus giving a meaningful expansion. Never- 
theless, by looking at the generic behaviour of these integral, one finds some evidence that 
it may converge. For instance, in the constant volume (M = 0) path, Eq. (3.29) gives 

A^^'=^-"(n;x„x^) = -^ f Al.e'^^kt+0,kt+e,k-+e,k-^^^^^^^-^j^^Q^^Q^-^ ^ (2n - l)/i(^3, ^4)) 



2TiGn J (2 



iTX 



(3.30) 

which is clearly finite (at least for n > 0). This is to be contrasted with the vacuum FRW 
case [2], and the example in [22], where a regulator is required in order to render finite the 
otherwise divergent amplitudes, even in this M = case. Further study of the vacuum case 
is in progress [23]. 
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IV. PROJECTION TO FRW 

At the end of section II, we commented on the projection from Bianchi I to FRW. We 
show here that when such projection is done at the level of the vertex expansion, Eq (3.18), 
one recovers the FRW vertex expansion of [2]. 

The structure in both cases is almost identical. One has a sum over volume sequences 
{um, ■ ■ ■ , Uq), and each amplitude is constructed by first obtaining a gravitational amplitude, 
and then performing the group averaging and scalar field integrations. Thus, all that remains 
is to show that the amplitude {xf\A{nM, • • • , "^o! ^m, ■ ■ ■ ,Ti', <^)\^i) given in Eq. (3.10) projects 
to the corresponding FRW one, 

A^^^iuM, . . . , no; TM, . . . , n; a) = e-^^^-^^^^r^M e^X-. >< 






• • • ^ ' ' ^nino 



when summing over all possible values of x/. 

To show this, it is convenient to explicitly write the intermediate anisotropics in the am- 
plitude {xf\A{nM, ■ ■ ■ ,no; tm, • • • , Ti; a) |xj). As before, this is done by introducing complete 
basis 1a = Xlx„ |x^)(x^| and 1a = Ey„ |ym)(ym| to the right and left of the A„^(r^+i-r™) 
operator in Eq. (3.10). Calling xj = jm and Xj = xq we have, 

^(yM|^(nM, •••, no; Tm,..., n; a)|xo) = ^ (yM|^nM(" - ^m)|xm) 

Ym xi,...,xm 

yoviYM 

{^M\VnMnM-i\yM~l){yM^l\AnM_^{TM -Tm-i)|xm-i) • • • (Xi |Kino IYo) (yol^m, (n) |xo) . 

(4.2) 

We now use the translation invariance of the operators to write each matrix elements in 
Eq. (4.2) as (x'|/|x) = (x' — x|/|0) where / is either the A or V operators. We then change 
the summation variables to y^ = y^ — x^ and x'^ = x^ — Jm-i- The different summations 
in Eq. (4.2) then decouple giving. 



^(yM|^(nM,---,no;rM,...,ri;a)|xo) = ^(y^^|A„^,(a - rM)|0) ^{^'mIVumum-iI^) 

^ ky'u-MnM-A^M - rM-l)|xM-l) ... ^(x'llKmolO) Xl^yol^no(^l)IO) 

,yW-i / \A J \y'o 

(4.3) 

Comparing Eq. (4.3) and Eq. (4.1) we see that our task reduces to showing that 
Ex(x|A„(r)|0) = e-^^®-"^ and J^M^WnlO) = 0^,^^ ■ That this is so can be seen as 
a direct consequence of the result in [19]. Let us nevertheless show it explicitly. 

For the V term, it suffices to look at Vn±in- The operator, given in Eq. (2.21) consist 
in an overall constant times six different translations in the x plane. Each term thus will 
pick up a single x from the sum. For instance, the first term gives a nonzero value only 
for X = (a^,0), in which case it gives a contribution of —ttG ^/n{n±T){2n ± 1). We then 
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conclude that 

$^(x|V;±i„|0) = -Q'nG^n{n±l){2n ± 1) = 0™ , (4.4) 

X 

as required. 

For the A„(r) term we have 

X 

Y, (k+|A+(r)|0)(k-|A;(r)|0) = 

k+,k-eZ2 

where in going from the first to second hne, we used Eq. (3.24). In going from the second 
to third hne, we used Eq. (3.25) and the identity X^fcez ^^'^^Z^'''" = <^(^) to directly evaluate 
the Fourier integrals. Using Eqs. (4.4) and (4.5) we then have 

^(x/|A(nM,---,no;rM,...,ri;a)|xi) = A^^^(nM, • • • ,no; tm, • • • ,ri; a) , (4.6) 

which implies 

y^ A{nM, • • • , no; Xj, x/; 0j, 0/) = A^'^'^{nM, • • • , no; 0i, 0/) . (4.7) 

Thus we see that the vertex expansion for our Bianchi model, Eq. (3.18) projects down to 
the vertex expansion of the FRW model, order by order. 

V. DISCUSSION 

Recently it has been shown [1] that one can take the Loop Quantum version of FRW 
cosmology and expand it as a sum over volume transition of amplitudes compatible with 
given initial and final states i.e. that the cosmological model of Loop Quantum Gravity can 
be re-written in terms of a sum over amplitudes, analogous to the spin foam approach. This 
sum over transition amplitudes is produced as a perturbation expansion of the constraint 
operator of LQC, thus linking perturbative dynamics of LQC to (the analogue of) spin 
foams. This analogue provided a useful link between the two theories, however because 
there is only one dynamic parameter in FRW cosmologies - the volume - the system has 
no analogue of the spin labels. In this paper we have extended the approach of [1] to 
the Bianchi I cosmological model, which, in addition to volume, has anisotropic degrees of 
freedom. We have shown that it is again possible to expand the dynamics of the model in 
terms of sums of amplitudes over volume transitions compatible with initial and final states. 
The additional anisotropic degrees of freedom of this model are analogous to the spin labels 
of spin networks, thus significantly improving the analogue to spin foams. 

The analogue remains at a formal level however, because one cannot directly associate the 
amplitudes with the changing of an underlying spin network. Despite this the association of 
the anisotropic degrees of freedom with the spin labels is well motivated by the fact that in 



14 

LQC they give the area (of the fiducial cell), which is precisely the role played by the spin 
labels (and edges) in a spin network. In addition to showing that the resulting summation 
over 'spin' labels is finite, we show that the projection to the FRW system occurs order by 
order in the expansion, thus recovering the results of [1]. 

Finally, although spin foams are typically taken to have spin changes only at vertices, 
it is generally expected that spin dynamics in the absence of graph changing vertices will 
play an important role in the final theory [24]. More precisely, that the action of the full 
constraint is non-trivial, even in the absence of vertices and hence that the amplitude for 
each vertex-free segment of the spin foam will be non-diagonal in the spin labels. In the 
analogue produced here we show that indeed the 'spin' changing amplitude is non-trivial, 
even in the absence of volume changing 'vertices'. Thus our full expansion is the analogue 
of a generalization of spin foams, allowing for 'spin' dynamics. 
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